Probing the Almeida-Thouless line away from the mean-field model 
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Results of Monte Carlo simulations of the one-dimensional long-range Ising spin glass with power- 
law interactions in the presence of a (random) field are presented. By tuning the exponent of the 
power-law interactions, we are able to scan the full range of possible behaviors from the infinite- 
range (Sherrington-Kirkpatrick) model to the short-range model. A finite-size scaling analysis of 
the correlation length indicates that the Almeida-Thouless line does not occur in the region with 
non-mean-field critical behavior in zero field. However, there is evidence that an Almeida-Thouless 
line does occur in the mean-field region. 
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I. INTRODUCTION 



The behavior of spin glasses in a magnetic field is 
still controversial. While the infinite-range (mean-field) 
Sherrington-Kirkpatrick (SK) modeli has a line of transi- 
tions at finite field known as the Almeida-Thouless (AT) 
line^ it has not been definitely established whether an 
AT line occurs in more realistic models with short range 
interactions. Previous numerical studiesi^SiLSiSiiSiii 
have yielded conflicting results: some data support the 
existence of an AT line in short-range spin glasses while 
others claim its absence. Recently^ a new approach 
using the correlation lcngthiSiiiiii^ at finite fields has 
been applied to the three-dimensional Edwards- Anderson 
Ising spin-glass model. The data of Ref. ^] indicate 
that, even for small fields, there is no AT line in three- 
dimensional spin glasses. 



Here we use the techniques of Ref. 11 to study the 
crossover between mean-field models and short-range 
spin glasses continuously by using a one-dimensional 
Ising chain with random power-law interactions. The 
model's advantages are twofold: first, a large range of 
system sizes can be studied, and second, by tuning the 
power-law exponent of the interactions, the universality 
class of the model can be changed continuously from the 
mean-field universality class to the short-range universal- 
ity class. We find that there appears to be no AT line 
for the range of the power-law exponent corresponding 
to a non-mean-field transition in zero field. By analogy, 
this suggests that there is also no AT line for short-range 
spin glasses, at least below the upper critical dimension. 



The paper is organized as follows. In Sec. [HJ we in- 
troduce in detail the model, observables, and numerical 
method used. In Sec. IIIII we present our results, and in 
Sec. IIVI we summarize our findings. 



II. MODEL, OBSERVABLES, AND 
NUMERICAL DETAILS 

The Hamiltonian of the one-dimensional long-range 
Ising spin glass with random power-law interactionsiSiil 
is given by 



n = - 



7 , JijSiSj 



hi Si 



(i) 



where Si = ±1 represents Ising spins evenly distributed 
on a ring of length L in order to ensure periodic boundary 
conditions. The sum is over all spins on the chain and 
the couplings J^- are given byi& 



(2) 



where the are chosen according to a Gaussian distri- 
bution with zero mean and standard deviation unity, and 
Tij = (L/tt) sin[(7r|i — j\)/L] represents the geometric dis- 
tance between the spins on the ring. 19 The power-law ex- 
ponent cr characterizes the interactions and, hence, deter- 
mines the universality class of the model. The constant 
c(o~) in Eq. @ is chosen to give a mean-field transition 
temperature T C MF = 1, where 



£[41 



c(a)' 



T — 



(3) 



Here [• • -] av denotes an average over disorder. 

In Eq. QJ, the spins couple to site-dependent random 
fields hi chosen from a Gaussian distribution with zero 
mean [/ij] av = and standard deviation [/if]iv 2 = #r- 
For a symmetric distribution of bonds, the sign of hi can 
be "gauged away" so a uniform field is completely equiv- 
alent to a bimodal distribution of fields with hi = ±Hr . 
While the AT line is usually studied for the case of a 
uniform field, the SK model with Gaussian random fields 
(as considered here) also shows an AT line. For short- 
range three-dimensional spin glasses it has been shown 
m Ref. [13 that results for Gaussian-distributed random 
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TABLE I: A summary of the behavior for different ranges 
of a in one space dimension and at zero field. 1R means 
infinite range, i.e., ^ii diverges unless the bonds Jij are 

scaled by an inverse power of the system size. LR means 
that the behavior is dominated by the long-range tail of the 
interactions, and SR means that the behavior is that of a 
short-range system. 

a behavior 



(7 = 

< a < 1/2 
1/2 < a < 2/3 
2/3 < a < 1 

1 < a < 2 
a > 2 



SK model 
1R 

LR (mean field with T c > 0) 
LR (non-mean field with T c > 0) 
LR (T c = 0) 

SR (T c = 0) 



fields agree within error bars with results for a uniform 
field. The use of Gaussian-distributed random fields has 
the advantage over a uniform external field in that we 
can apply a useful equilibration test?ii*i2i2£ see Eq. <|1(J[) 
below. 

From Eq. ®, we see that for a < 1/2, c(a) varies 
with a power of the system, c(cr) ~ £-( 1 - 2cr )/ 2 ; for large 
L. We shall denote systems in this region as "infinite 
range" (IR). The extreme limit of a = gives the SK 
model, whose solution is the mean-field (MF) theory for 
spin glasses. For a > 1/2, c(a) tends to a constant as 
L — > oo. As discussed in an earlier work (see Ref . Il8l and 
references therein), for 1/2 < a < 1, the system has a 
finite-temperature transition into a spin-glass phase in a 
long-range (LR) universality class at zero field. For 1 < 
a < 2, the system has T c = and the critical behavior is 
also determined by the LR universality class. For a > 2, 
we have a short-range (SR) universality class with T c = 0. 
Finally, we note 21 that for 1/2 < a < 2/3 the critical 
behavior is mean- field-like, while for 2/3 < a < 1 it 
is non-mean field like. This behavior is summarized in 
Table [I] Critical exponents depend continuously on a 
in the LR regime, but are independent of er in the SR 
region. Here we focus on the regime 1/2 < a < 1 because 
there the system exhibits a finite-temperature transition 
that can be tuned continuously away from the mean-field 
universality limit by changing the exponent a. 

To determine the existence of an AT line, we compute 
the two-point correlation Iengthiii*i2ii4 We calculate the 
wave-vector-dependent spin-glass susceptibility which is 
denned by 



1 



Xsc(k) = -^ 



{SiSj)T — (Si)T(Sj)T 



ik-(Ri-Rj) 



(4) 

where (■ ■ -)t denotes a thermal average. Note that at 
zero field (Si)x can be set to zero. The correlation length 
of the finite system is then given by 



1 



2 sin(/c min /2) 



Xsg(O) 

XSG(k m in 



- 1 



1/(2,7-1) 



(5) 



-0.5 



mr| — i i | — i i i i — i i mi — i i inn — i i i i — r 
mj ooooooooo oo o_o_8-a- a a a 



1.5 



-2.5 




10 a 10 3 10 4 10 5 10 6 10 7 



FIG. 1: (Color online) Sample equilibration plot for a — 0.85, 
L = 128, _Hr = 0.1, and T — 0.10 (the lowest temperature 
simulated at finite fields). Data for the average energy U, 
and U(qi, q) defined in Eq. 1101 . as a function of equilibration 
time i eq - They approach their common value from opposite 
directions and, once they agree, do not change on further 
increasing t cq . The inset shows data for the correlation length 
divided by system size as a function of equilibration time. The 
data are independent of t e q once U and U (qi , q) agree. 



where k m ; n = (27r/L,0, 0) is the smallest nonzero wave 
vector. The reason for the power l/(2cr — 1) is that 
at long wavelengths, we expect a modified Ornstein- 
Zernicke forr 



Xsc(k) oc (v + k 



(G) 



for the long-range case, where v is a measure of the devi- 
ation from criticality. It follows that the bulk correlation 
length £ diverges for v — > like v^ 1 ^ 2rT ~ 1 \ 

The correlation length divided by the system size £,l/L 
has the following scaling property: 



^- = x(l 1 / v [t-t c {Hk)] 



(7) 



where v is the correlation length exponent and T c (Hn) 
is the transition temperature for a field of strength Hr . 
This behavior is similar to that of the Binder ratio^ but 
it shows a clearer signature of the transition as the data 
are not restricted to a finite interval. 

In order to test equilibration of the Monte Carlo 
method, we also compute the link overlap^ qi given by 



qi ~ ]v ^ (T c Mi? ) 2 [{ 1 j 1 j ' 



Tjav , 



(8) 



where T C MF is given by Eq. © and a and (3 refer to 
two replicas of the system with the same disorder. In 
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TABLE II: Parameters of the simulations for — 0.0. 
A sa is the number of samples, A sw is the total number of 
Monte Carlo sweeps for each of the 2At replicas for a single 
sample, T m - m is the lowest temperature simulated, and Nt is 
the number of temperatures used in the parallel tempering 
method for each system size L and power-law exponent a. 



(7 


L 


A sa 


A sw 


Train 


N T 


0.55 


32 


5000 


10240 


0.405 


15 


0.55 


64 


5000 


10240 


0.405 


15 


0.55 


128 


5000 


20480 


0.405 


15 


0.55 


256 


5000 


102400 


0.405 


15 


0.55 


512 


5000 


32768 


0.630 


11 


0.65 


32 


5000 


10240 


0.405 


15 


0.65 


64 


5000 


10240 


0.405 


15 


0.65 


128 


5000 


20480 


0.405 


15 


0.65 


256 


5000 


102400 


0.405 


15 


0.65 


512 


5000 


524288 


0.405 


15 


0.75 


32 


5000 


10240 


0.405 


15 


0.75 


64 


5000 


10240 


0.405 


15 


0.75 


128 


5000 


20480 


0.405 


15 


0.75 


256 


5000 


102400 


0.405 


15 


0. to 


512 


2500 


524288 


A A AC 
0.4UO 


lo 


0.85 


32 


5000 


10240 


0.405 


15 


U.oO 


64 


5000 


20480 


U.4Uo 


1 K 
1Q 


0.85 


128 


5000 


102400 


0.405 


15 


0.85 


256 


5000 


204800 


0.405 


15 


0.85 


512 


2500 


204800 


0.405 


15 


addition, 


we compute the spin overlap q 


given by 








1 

q= N 


N 

^[(SfSfHav. 

i=l 




(9) 



Because both the fields and interactions have a Gaus- 
sian distribution, integrating by parts the expression for 
the average energy per spin U givesiiiiSiSSi 



U=U(q u q) 



(TP 2 
2T 



Hi 
T 



(1-?). (10) 



As shown in Fig.^ when starting from a random spin 
configuration, U approaches its equilibrium value from 
above while U (qi , q) approaches its equilibrium value 
from below. Once U — U(qi,q), the data do not change 
by further increasing the number of Monte Carlo steps, 
which shows that the system is in equilibrium. It is also 
important to ensure that other observables are also in 
equilibrium once U = U(qi,q), and this is shown in the 
inset to Fig. ^ for the case of the correlation length. 

The simulations are done using the parallel tempering 
Monte Carlo methodiS&Si The method is not as efficient 
in a field, 25,26 but nevertheless it performs considerably 
better than simple Monte Carlo. In order to compute 



TABLE III: Parameters of the simulations for J/r = 0.1. 
N sa is the number of samples, A sw is the total number of 
Monte Carlo sweeps for each of the 4At replicas for a single 
sample, T m i n is the lowest temperature simulated, and At is 
the number of temperatures used in the parallel tempering 
method for each system size L and power-law exponent a. 



a 


L 


A sa 


A sw 


Tmin 


A T 


0.55 


32 


5000 


81920 


0.100 


26 


0.55 


64 


5000 


327680 


0.100 


26 


0.55 


128 


5000 


1310720 


0.100 


26 


0.55 


256 


2000 


1048576 


0.405 


15 


0.55 


512 


2000 


65536 


0.760 


9 


0.65 


32 


5000 


81920 


0.100 


26 


0.65 


64 


5000 


327680 


0.100 


26 


0.65 


128 


5000 


1310720 


0.100 


26 


0.65 


256 


2000 


1048576 


0.195 


20 


0.65 


512 


2000 


524288 


0.500 


13 


0.75 


32 


5000 


81920 


0.100 


26 


0.75 


64 


5000 


327680 


0.100 


26 


0.75 


128 


5000 


1310720 


0.100 


26 


0.75 


256 


2000 


8388608 


0.100 


26 


0.85 


32 


5000 


81920 


0.100 


26 


0.85 


64 


5000 


327680 


0.100 


26 


0.85 


128 


2000 


16777216 


0.100 


26 



the products of up to four thermal averages in Eq. Q 
without bias, we simulate four copies (replicas) of the 
system with the same bonds and fields at each temper- 
ature. Simulations are performed at zero field, as well 
as at = 0.1, a field that is considerably smaller than 
the er-dependent transition temperature T c . Parameters 
of the simulations at zero and finite fields are presented 
in Tables ITTI and Hill respectively. 



III. RESULTS 

We first consider the case of zero field and take a = 
0.55, 0.65, 0.75, and 0.85. The values a = 0.75 and 
0.85 are in the non-MF region (see Table [IJ while a = 
0.55 is in the MF region and, furthermore, is close to 
the value (a — 1/2) where the system becomes infinite 
range. The value a = 0.65 is close to the point a = 2/3 
where the critical behavior changes from MF to non-MF. 
The data are shown in Fig. [5] In all cases, the data 
cross at a transition temperature which we determine as 
T c = 1.03(3) for cr = 0.55, 0.86(2) for a = 0.65, 0.69(1) 
for a = 0.75, and 0.49(1) for a = 0.85. Note that T c 
decreases continuously with increasing a and is expected 
to drop to zero at a = 1. 18 For the SK model (cr = 0), one 
has T c = 1, essentially the result we find for a = 0.55, 
so it is possible that T c has little variation with a for 
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T T 

FIG. 2: (Color online) Each figure shows data for (,l/L vs T for Hr = for different system sizes, for a particular value of a. 
For all values of a, the data cross indicating that there is a spin-glass transition at finite temperature. 



a < 0.55. 

Next we consider = 0.10 and show the data in 
Fig. The results for a = 0.75 and 0.85, which are in 
the non-mean-field regime, show no sign of a transition. 
However, the data for a = 0.55 do show a signature of a 
transition at T c = 0.96(2) < T C (H R = 0). Whether this 
would persist up to infinite system sizes is not clear, but 
it certainly cannot be ruled out. The results for a — 0.55 
show that the method used here is capable of detecting 
an AT line in the presence of a field. For a = 0.65, 
the data shows a marginal behavior. Since a = 0.65 is 
close to the value of 2/3 which separates MF and non- 
MF behavior in zero field, this marginal behavior may 
indicate that 2/3 is also the borderline value below which 
an AT line occurs. An alternative possibility, which we 
cannot rule out, is that an AT line only occurs in the 



infinite-range region (a < 1/2) but that as a is decreased 
toward 1/2, one needs to study larger system sizes to sec 
the absence of a transition. 



IV. CONCLUSIONS 

We have considered a one-dimensional spin-glass model 
with long-range interactions that allows the universality 
class to be changed from the infinite-range limit to the 
short-range case by tuning the power-law exponent a of 
the interactions. We find that there does not appear to 
be an AT line in a field for models with a in the range 
where there is non-mean-field critical behavior at zero 
field. However, in the region of a that is not infinite-range 
but has mean-field critical behavior (1/2 < a < 2/3), 
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T T 

FIG. 3: (Color online) Each figure shows results for (i/LvsT for Hb. = 0.10 for different system sizes, for a particular value of 
a. For a = 0.55, the data cross at T c = 0.96(2) showing that an AT line seems to be present. For a — 0.65, the data show close 
to marginal behavior. This may indicate that a = 0.65 is close to the borderline value for having an AT line. For a — 0.75 
and a = 0.85, the data do not cross for any temperature down to T — 0.10, which is considerably smaller than the zero field 
transition temperatures. This indicates that there is no AT line. Overall the results indicate that, on increasing a, the AT line 
disappears. 



there does appear to be an AT line. 

These conclusions rely on extrapolating from finite 
sizes to the thermodynamic limit. It would be partic- 
ularly interesting to know if the conclusion that there is 
an AT line for a = 0.55 persists in the thermodynamic 
limit, or whether the AT line really only occurs in the 
infinite-range case (er < 1/2). It is possible that as a is 
decreased, larger sizes or larger values of £/r are needed 
to probe the asymptotic behavior. Therefore, it would be 
desirable to simulate a range of values of Hr, especially 
for a = 0.55. We have some results for £/r = 0.2 for 
a relatively small number of samples and for sizes only 



up to L = 96, which indicate a crossing at a lower tem- 
perature than for H-r = 0.1. However, we are unable 
to carry out a systematic study of the dependence on 
.Hr because the results presented above already required 
considerable computer time, and the parallel tempering 
algorithm becomes less efficient at larger fields. 

Making an analogy between the one-dimensional long- 
range model for different values of a and short-range 
models for different values of space dimension d, we infer 
that there is no AT line for short-range spin glasses in 
the non-mean-field regime, i.e., below the upper critical 
dimension d u = 6. However, there may be an AT line 
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above the upper critical dimension. Speculations along 
these lines have also been made very recently by Moore p2L 
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